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I. INTRODUCTION 



The state space of an n-dimensional quantum system is the set of aU nx n positive semidefinite complex matrices of 
trace 1 cahed density matrices. The density matrices of rank one (pure states) can be identified with nonzero vectors 
in a complex Hilbert space of dimension n. However, one has to take into account that the same state is described 
by a vector ip and Xip, where A 0. Hence pure states are in one-to-one correspondence with rays {A^ : ^ X € C}. 
The rays form a smooth manifold called complex projective space CP"~^. In the infinite-dimensional case we have 
to consider density operators instead of density matrices, and the space of pure states is the complex projective space 
over the infinite-dimensional Hilbert space. 

The problem of measuring a distance between two quantum states with a suitable metric attracts a lot of attention 
in recent years. The Hilbert-Schmidt norm of an operator ||A||2 = ■\/Tr(At A) induces a natural distance between 

two density operators dnsipi, P2) = v^Tr [(pi — p2)'^]- This distance has been recently used to describe the dynamics 
of the field in Jaynes-Cummings model M and to characterize the distance between certain states used in quantum 
optics Another distance generated by the trace norm ||^||i — TrV A^A was used by Hillery to measure the 
non-classical properties of quantum states. 

A concept of statistical distance was introduced by Wootters in the context of measurements which optimally 
resolve neighboring pure quantum states. This distance, leading to the Fubini-Study metric in the complex projective 
space, was later generalized to density matrices by Braustein and Caves Q and its dynamics for a two-state system 
was studied by Braunstein and Milburn 0. It was shown |^ that for neighboring density matrices the statistical 
distance is proportional to the distance introduced by Bures in late sixties The latter was studied by Uhlmann 
and Hiibner [|o|, who found an explicit formula for the Bures distance between two density operators P2) = 

2(l — iT:[{p\^'^P2p\^'^Y^'^\). Note that various Riemannian metrics on the spaces of quantum states were also considered 
by many other authors (see |pT|~p^). 

In the present paper we introduce yet another distance in the space of density operators, which fulfills a following 
semiclassical property : the distance between two coherent states |ai) and \a2) localized at points ai and 02 of the 
classical phase space n endowed with a metric d is equal to the distance of these points 

D{\ai),\a2))^d{aua2). (1.1) 

This condition is quite natural in the semiclassical regime, where the quasi-probability distribution of a quantum 
state tends to be strongly localized in the vicinity of the corresponding classical point. A motivation to study such 
a distance stems from the search for quantum Lyapunov exponent, where a link between distances in Hilbert space 
and in the classical phase space is required [ ^ . 

In order to find a metric satisfying conditio n (|l.l[ ) it is convenient to represent a quantum state p in the Q- 
representation (called also the Husimi function) ||l^ 

Hp{a) {a\p\a), (1.2) 

defined with the help of a family of coherent states |a) , a G fi, which fulfill the identity resolution 1 = \a) {a\dm{a), 
where m is the natural measure on fi. For a pure state lip) one has H^{a) = |(V'|Q!)p. The choice of coherent states 
is somewhat arbitrary and in fact one may work with different systems of coherent states . In the present paper 
we use the classical harmonic oscillator (field) coherent states, where ft = C and dm{a) — d'^a/n. For convenience 
we will use in the sequel the renormalized version of the Husimi function putting H{a) = {a\p\a)/TT and integrating 
over d'^a. Note that the Husimi representation of a given density operator determines its uniquely [^. Since Husimi 
distributions are non-negative and normalized, i.e., J^H{a)dm{a) = 1 and _ff > 0, it follows that a metric in the 
space of probability densities Q : 17 — > i?+ induces a metric in the state space. 

In this work we propose to measure the distance between density operators by the Monge distance between the 
corresponding Husimi distributions. The original Monge problem consists in finding an optimal way of transforming 
a pile of sand into a new location. The Monge distance between two piles is given by the 'path' traveled by all grains 



under the optimal transformation |E0 21 



This paper is organized as follows. In sect. II we give a definition of the Monge metric, present an explicit solution for 
the one-dimensional case and discuss some methods of tackling the two-dimensional problem. Sect. Ill contains some 
examples of computing the Monge distance between certain states often encountered in quantum optics. Concluding 
remarks are given in sect. IV. Variational approach to the Monge problem is briefly presented in the Appendix. 



II. MONGE METRIC 
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A. Monge transport problem 



The original Monge problem, formulated in 1781 |22j, emerged from studying the most efficient way of transporting 
soil [|0|: 

Split two equally large volumes into infinitely small particles and then associate them with each other so that 
the sum of products of these paths of the particles to the volume is least. Along what paths must the particles be 
transported and what is the smallest transportation cost? 

Fig. 1 represents a scheme for this problem. Here we give a general definition of the Monge distance between two 
probability densities Qi and Q2 defined in 5* = i?". Let fii and SI2, determined by Qi, describe the initial and the 
final location of 'sand': = {{x,y) E S x _R+ : < y < Qi{x)}. Due to normalization of Qi the integral d^xdy 
is equal to unity. Consider maps T : S S which generate volume preserving transformations ili into i.e., 

QUx)^Q2iTx)\r{x)\ (2.1) 

for all X G S, where T'{x) denotes the Jacobian of the map T at point x. We shall look for a transformation giving 
the minimal displacement integral and define the Monge distance pO| , pT[ | 

Dm{Qi,Q2) ■■= inf [ \x- T{x)\Qi (x) d^'x , (2.2) 
Js 

where the infimum is taken over all T as above. The optimal transformation (if exists) Tm is called a Monge plan. 
Note that in this formulation of the problem the "vertical" component of the sand movement is neglected. 

It is easy to see that Monge distance fulfills all the axioms of a metric. This allows us to define a 'classical' distance 
between two quantum states pi and p2 via the Monge distance between the corresponding Husimi distributions: 

Dci{puP2) ■.= DM{Hi{a),H2{a)). (2.3) 

The Monge distance satisfies the semiclassical property: it is shown below that the distance between two coherent 
states, represented by Gaussian Husimi distributions localized at points oi and 02, equals to the classical distance 
|ai - a2\. 

It is sometimes useful to generalize the notion of the Monge metric and to define a family of distances -Djv/p labeled 
by a continuous index p (0 < p < 00) in an analogous way: 

[Dm,{Qi,Q2)Y ■=mi I \x-T{x)\PQi{x)d^x . (2.4) 
Js 

For p — 1 one recovers the Monge distance I?a/i — Dm, while the Frechet distance is obtained for p 2. A 
more general approach to the transport problem was proposed by Kantorovitch and further developed by Sudakov 
[Q. In contrast with the definition of Monge discussed in this work, the Kantorovitch distance between Qi and Q2 
is explicitly symmetric with respect to exchange of both distributions. For a comprehensive review of metrics in the 
space of probability measures and other generalizations of the Monge distance see the monograph of Rachev . 

B. Salvemini solution for ID problem 

For S = R the Monge distance can be expressed explicitly with the help of distribution functions Fi{x) — 
Qi{t)dt, i — 1,2. Salvemini |^ and Dall'Aglio |2^ obtained the following solution of the problem for p = 1: 

/+00 
\F,ix) - F2ix)\dx , (2.5) 
-00 

represented schematically in Fig. 2. 

This result was generalized in the fifties to all p > 1 by several authors (see pQ,pil ) . We have 



[Dmp{Qi,Q2W = / \Fr\t) - F2\t)\Pdt . (2.6) 
Jo 



3 



C. Poisson-Ampere-Monge equation for 2D case 



Consider smooth densities Qi, Q2 ■ ^ R^- We are looking for a transformation field w ~ (^1,^2) : R^ ^ R^ 
fulfilling wlxi, X2) — Tm{xi, X2) — {xi, X2), where Tm is an optimal Monge transformation minimizing the right-hand 
side of (|]^) • Restricting ourselves to smooth transformations T we may apply the standard variational search for the 
optimal field w. In Appendix it is shown that vot{w) = if p = 2, and the potential ip : w = grad((/j) satisfies the 
following partial differential equation 

_ ( )2 _ Ql{xi,X2) , s 

Q2(Xl + (Pxi,X2 + •^X2) 



Solving this Laplace-Ampere-Monge (LAM) equation for the potential ip we get the Frechet distance Dm2 from (2.4) 
computing the minimal displacement 

/-t-00 r+00 
dxi / dX2{^l,+^ljQl{xi,X2) . (2.8) 
-00 — 00 



Though it is hardly possible to solve equation (2.7) in a general case, it might be used to check, whether a given 



transformation can be a solution of th e M onge problem. Let us remark that an optimal transformation field w (if 



exists) must not be unique. Moreover, (2/7) provides only a sufficient condition for w being optimal. It is important 
to note that the two dimensional Monge problem posed two hundreds years ago has not been solved in a general case 
until now iol,^. 



D. Estimation of the Monge distance via transport problem 



One of major tasks of linear programming is the optimization of the following transport problem. Consider N sup- 
pliers producing ai {i = 1, . . . , N) pieces of a product per a time period and M customers requiring bj {j — I, . . . , M) 
pieces of the product at the same time. Let [cij) he a, N x M cost matrix, representing for example the distances 
between sites. Find the optimal transporting scheme, minimizing the total transport costs C = X^il^i X^^li ^ij^ij- 
The non-negative elements of the unknown matrix (xy ) denote the number of products moved from i-th supplier to 
j-th customer. The optimization problem is subjected to the following constrains: X^^li •'^ij ~ ^^'^ ^f=i ^ij ^j- 

In the simplest case the total supply equals the total demand and ~ S^li ^j- 

The transport problem described above gave, with a related assignment problem, an impulse to develop methods 
of linear programming already half a century ago |2^ , |2^ . Since then several algorithms for solving the transport 
pro blem numerically have been proposed. Some of them, as the northwest corner procedure and Vogel's approximation 
p8| , are available in specialized software packages. It is worth to add that the transport methods are widely used to 
solve a variety of problems in business and economy such as, for instance, market distribution, production planning, 
plant location and scheduling problems. 

It is easy to see that the transport problem is a discretized version of the Monge problem defined for continuous 
distribution functions. One can, therefore, approximate the Monge distance between two distributions Qi{x) and 
Q2ix) (where x stands for a two dimensional vector), by solving the transport problem for delta peaks approximation 
of the continuous distributions: qi = '^j^^iQi{xi)S{x — Xi) and (72 — '^j=iQ2ixj)6{x — Xj). The quality of this 
approximation depends on the numbers TV and M of peaks representing the initial and the final distribution, respec- 
tively, and also on their location with respect to the shape of both distributions. Numerical study performed with 
the nortiiwest corner algorithm for some analytically soluble examples of the 2-D Monge problem shows that for 
reasonably smooth distributions one obtains Monge distance with fair accuracy for a number of peaks of the order of 
hundreds. 



III. MONGE DISTANCE BETWEEN SOME STATES OF QUANTUM OPTICS 

In this section we compute Monge distances between certain quantum states. Even though our considerations are 
valid in the general framework of quantum mechanics, for the sake of concreteness we will use the language of quantum 
optics. Let us recall briefly the necessary definitions and properties. 

Let a and be the annihilation and creation operators satisfying the commutation relation [a, a^] — 1. The 'vacuum' 
state |0) is distinguished by the relation a|0) = 0. Standard harmonic-oscillator coherent states \a) can be defined as 
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the eigenstates of the annihilation operator a\a) = a\a) or by the Glauber translation operator D{a) = exp{aa^ — a* a) 
as \a) = D{a)\0). Coherent states, determined by an arbitrary complex number a = Xi -\- 1x2, enjoy a minimum 
uncertainty property. They are not orthogonal and do overlap. The Husimi distribution of a coherent state |/3) is 
Gaussian 

Hp{a) = -|(/3|a)P = ^ exp(-|a - /Jp). (3.1) 

TT TT 

Squeezed states \j,a) also minimize the uncertainty relation, however the variances of both canonically coupled 
variables are not equal. They are defined as 

\^,a):^D{a)S{^)\Q), (3.2) 

where the squeezing operator is 5'(7) — exp[i(7*a^ — ja^'^)]. The modulus g of the complex number 7 = ge^^^ 
determines the strength of squeezing, s — — I, while the angle 9 orients the squeezing axis. The Husimi distribution 
of a squeezed state \"f,P) is a non-symmetric Gaussian localized at point (3 and for 9 = reads 

H^A^i^X2) = - exp[-(Re(/3) - xif/{.s + if - (Im(/3) - .T2)'(.s + 1)^] (3.3) 

TT 

Each pure state can be expressed in the Fock basis consisting of n-photon states \n), n = 0, 1, 2, . . .. Each coherent 
state can be expanded in the Fock basis as 

a 



= p-i"r/2 



E^l")- (3-4) 



n=0 

The known scalar product {a\n) allows one to write the Husimi function of a Foci state 

i?|n)(«) = -^e-H . (3.5) 

The vacuum state |0) can be thus regarded as the single Fock state being simultaneously coherent. 

In contrast to the above mentioned pure states, the thermal mixture of states with the mean number of photon 
equal to n is represented by the density operator 



'"^I]?^— TvI+tI")("I- (3-6) 



n=0 

Its Husimi distribution is the Gaussian centered at zero with the width depending on the mean photon number. 



A. Two coherent states 



Let us consider an arbitrary two dimensional distribution function Qi{xi, X2) and a shifted one Q2{xi, X2) = 
Qiixi — xl, X2 — X2)- It is intuitive to expect that the simple translation given by w{xi, X2) = const = (a;J, X2) solves 
the corresponding Monge problem. Since for the respective potential (j){xi, X2) = x\xi + X2X2 the second derivatives 
vanish, the both two sides of the LAM equation (2.7) are equal to zero and the maximization condition is fulfilled. 

It follows from this observation that for two coherent states defined on the complex plane the Monge plan reduces 
to translation. Integration in ( ^.2| ) is trivial and provides the Monge distance between two arbitrary coherent states 
I a) and 

D,i{\a),m = \a-(i\. (3.8) 

This is exactly the semiclassical property we demanded from the metric in the state space. The distance of the 
coherent state |a) from the vacuum state |0) is equal to \a\, which simply is the square root of the mean number of 
photons in the state \a). The classical property is fulfilled by the generalized Monge metric -Djv/p for any positive 
parameter p. 
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B. Coherent and squeezed states 



We compute the distance between a coherent state |a) and the corresponding state squeezed I7, a). Due to invariance 
of the Monge optimization with respect to translations this distance is equal to the distance between vacuum |0) and 
the squeezed vacuum |7,0). For simplicity we will assume that squeezing is performed along the xi axis, i.e., the 
complex squeezing parameter is real ^ = g E R. 

The corresponding Monge problem consists in finding the optimal transformation of the symmetric Gaussian 
Qi{xi,X2) = exp(— a;^ — a^l)/"' i^^to an asymmetric one Q2ixi,X2) = exp(— a;^/(s + 1)^ — 2:2(3 + l)^)/7r. Consid- 
ering contours of the Husimi distribution, often used to represent a state in quantum optics, a circle has to be 
transformed into an ellipse. If p = 2 then the following affine transformation T{xi,X2) = {xi/{s + 1), 2:2(5 + 1)) 
corresponds to the irrotational transport field w(xi, 22) = {—sxi/{s + 1), S22). It can be obt aine d as the gradient of 



the potential (p{xi, 22) = —sxi/{2s + 2) + SX2/2, for which both sides of the LAM equation (2.7) vanish. Hence field 
w provides a Monge plan for this problem and the distance is given by its length integrated over the volume of Qi 

D,m, \g, 0» = - / d2i / dx2 exp(-2? - xl)Jxll{s + 1)2 +xl^ —E , (3.9) 

where E{x) stands for the complete elliptic integral of the second kind. 
A simpler result may be obtained in this case for the Frechet distance 

[i?M.(|0),|.g,0))] / d2i / dx2 exp(-2:^ - xl){xl/{s + 1)^ + x^) = -(1 + ). (3.10) 

Eventually, the 'Chebyshev' distance I^a/^, characterizing the maximal dislocation length, is equal to s. For large 
squeezing all three distances grow proportionally to the squeezing parameter s, with slopes k\ = ka = l/v^ 

and kao = 1 ordered according to the index p. This agrees with an observation that pi < P2 implies Dm < Dm 

13- 

C. Vacuum and thermal states 

Since Husimi distributions of both states is rotationally invariant, it is convenient to use radial components of 
the distribution Ri{r) = 27rr_ffi(r, </)). One may then reduce the problem to one dimens ion and find the Monge 



distance via radial distribution functions Fi{r) — Ri{r')dr' . Taking Husimi distributions (3T) and (3/7) we get the 



corresponding distribution functions Fi{r) = 1 — cxp(— r2) and F2{r) = 1 — exp(— r2/(rt + 1)). Using the Salvcmini 



formula (2.5) we obtain 

/>C30 / 

D,m{OlPn)^J^ \F,{r)-F2{r)\dr^^{V^-l). (3.11) 
In the same way we get a more general formula for the Monge distance between two thermal states 



F>ci{Pni,Pn2) = —\Vni + 1 - V"-2 + l|- (3-12) 



D. Two Fock states 

A s in the previous example the rotational symmetry of the Fock states allows us to use the ID formula. Integrating 
( |3.5| ) for a Fock state \n) we can express the radial distribution function in terms of the incomplete Gamma function 
r(2:, r) as 

r(n + l,r2) 

^"(") = ^- r(n + i) ■ (^-'^^ 



Since for different Fock states the distribution functions do not cross, applying Salvemini formula (2^) we obtain the 
Monge distance 



D,i{\m),\n)) = 



pOO pGO 

\ Fn{r)dr- F.rn{r)dr\ = V^\Crn - Cn\, (3.14) 
Jo Jo 



where the integrals Ci can be found analytically: Co = 1/2, Ci = 3/4, C2 15/16, C3 = 35/32, . 
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IV. DISCUSSION 



We have presented a definition of distance between quantum states (i.e. elements of a Hilbert space) which possesses 
a certain classical property, natural for investigation of the semiclassical limit of quantum mechanics. Monge distance 
between the corresponding Husimi functions fulfills the axioms of a metric and induces a 'classical' topology in the 
Hilbert space. It is worth to emphasize that the Monge distance between two given density matrices depends on the 
topology of the corresponding classical phase space. 

Consider a quantum state prepared as a superposition of two coherent states separated in the phase space by x. It is 
known [^o|-^ that such a state interacting with an environment evolves quickly toward a mixture of the two localized 
states. The decoherence time decreases with the classical distance x, equal just to the Monge distance between both 
coherent states. We expect therefore that the Monge distance between two arbitrary quantum states might be useful 
to determine the rate with which the superposition of these two states suffers the decoherence. 

It is possible to generalize our approach in several directions. Instead of the standard Husimi distributions used 
throughout this paper, one may study Monge distances between generalized Husimi distributions Hp{a) — {a\p\a) 



where \a) are generalized coherent states [ p3| . For example, one may use for this purpose squeezed states |34 , or the 
spin coherent states jH,^, if the classical space is the two-dimensional sphere. 

Moreover, our considerations based on the Husimi representation of quantum states, may be in fact extended to the 
Wigner function. In spite of the fact that the Wigner function can take on negative values, it is normalized and the 
Monge problem of finding an optimal way to transport one Wigner function into another might also be considered. 
The concept of the classical distance between two Husimi (Wigner) functions is not only of theoretical interest, since 
novel methods of measuring Husimi and Wigner distributions have been recently developed 
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APPENDIX A: VARIATIONAL APPROACH TO MONGE PROBLEM 



1. ID case 



Let Qi and Q2 be smooth densities. In ID case there is only one map T fulfilling ( p.l[ ). It can be described with the 
aid of distribution functions as T{x) = F2^{Fi{x)), where Fi{x) = Qi{t)dt for x £ R. This allows us to express 
the generalized Monge distance as an integral (|2^ 



[Dmp{Qi,Q2)V 



Qiix)\F2^[Fi{x)]~x\Pdx. 



(Al) 



In the simplest case p = 1 (Al) reduces to the Salvemini formula (|2.5|) and for p > 1 to formula (2.6) 



2. 2D case 



Consider two smooth densities Qi(a;i, X2) and Q2ix i, X 2)- We are looking for a map T{xi, X2) = (a;-|-wi(a;i, X2), X2 + 
W2 (2^1,2^2)) transforming Qi into Q2 (i.e. such that ( ^.l[ ) is fulfilled) and minimizing the quantity 

Ql{xi,X2)\wl{xi,X2) + 'W2ixi,X2)\''^^dxidX2 ■ (A2) 

-00 

The index p, labeling the generalized distance, is equal to one for the Monge metric and to two for the Frcchct metric. 
Introducing a Lagrange factor A we write the Lagrange function in the form 

Lp - Qiiwl + wlfl^ + A(Qi - 02(T)) [(1 + wi, J(l + W2^^) - w^^^W2.^\ . (A3) 

The Lagrange-Euler equations for variations of Lp allow us to obtain the partial derivatives of A 
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Aj;, = 2pCp{wi{l + wix^) + W2W2X1) 
X^^ = 2pCp(w2(l + W2X2) + W1W1X2) 



■^xiX2 = ^X2xi we get 



where Cp = {wf + W2YP ^^Z^. Using the equahty 

Wix2 {wlip- 1) + W2) - W2xt {wlip - 1) + wj) + ip- 2){W2X2 ~ WixJwiW2 = 



(A4) 



(A5) 



li p — 2 we deduce from ( [A5|) that wix2 — W2xi, i-G-j rot(i/;) = 0. Taking the potential ip : w — grad(93) we see that 
ip fulfills (2/7) and formula (|2.8D holds. 
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FIG. 1. Monge transport problem: How to shovel a pile of sand Qi{xi,X2) into a new location Q2{xi,X2) minimizing the 
work done? 



FIG. 2. Monge distance between ID functions Qi{x) and Q2ix) may be represented as the area between graphs of the 
corresponding distribution functions Fi (x) and F2 (x) 
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